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^>«^ ' We present an analytical and numerical study of the Klein-Gordon kink-soliton dynamics in 

inhomogeneous media. In particular, we study an external field that is almost constant for the whole 
system but that changes its sign at the center of coordinates and a localized impurity with finite- 
Pi^' width. The soliton solution of the Klein-Gordon-like equations is usually treated as a structureless 
, point-like particle. A richer dynamics is unveiled when the extended character of the soliton is taken 
into account. We show that interesting spatiotemporal phenomena appear when the structure of the 
soliton interacts with finite-width inhomogeneities. We solve an inverse problem in order to have 
external perturbations which are generic and topologically equivalent to well-known bifurcation 
models and such that the stability problem can be solved exactly. We also show the different 
quasiperiodic and chaotic motions the soliton undergoes as a time-dependent force pumps energy 
into the traslational mode of the kink and relate these dynamics with the excitation of the shape 
modes of the soliton. 
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I. INTRODUCTION 



j— * ■ For a variety of systems the interplay between nonlinearity and disorder results in novel and fascinating phenomena 
" 10,^ . Particularly, the study of soliton dynamics in inhomogeneous and disordered media has received a great deal of 
Oh' attention in recent years since it concerns real condensed matter systems and phenomena. 

^ ' It is well known that transport properties in inhomogeneous and disordered media can change dramatically when the 
nonlinearity allows the creation of solitons. As a first step in the study of soliton dynamics in disordered media many 
authors explored the interaction between a soliton and an isolated impurity. This soliton-impurity interaction has been 
modeled by a point-like impurity and a structureless soliton. In real situations we can have several inhomogeneities of 
different kinds. When the distance between impurities is considerable higher than the width of solitons and impurities 
the traditional approach gives a correct result. When this is not so, we witness a series of surprising phenomena. 

The dimensionless Klein-Gordon- like equations model a wide variety of soliton bearing systems 0J|,|ll|,|l^-0 , 
including charge density waves, Josephson junctions, structural phase transitions, crystal growth, polymers, proton 
conductivity, macromolecules and hydrogen-bond chains: 

(t)xx - 't'tt - ^ ^ J^{4>,4>t,x,t); (1) 

here U — U{(j)) is a potential that possesses at least two minima [Q, meanwhile T{(t), (pt, x,t) represent additional 
forces (external forcing, dissipation, presence of impurities, inhomogeneous external fields, coupling to other degrees 
of freedom) . 

The soliton solution of Eq. (^ is usually treated as a structureless point-like particle. A richer dynamics is unveiled 
when the extended character of the soliton is taken into account. For instance, in a previous work |^], Gonzalez and 
Holyst studied the (j)'^ equation {U = |(0'^ — 1)^): 

- 0« - T'/'t + ^ (0 - = -F{x) - G{x, t). (2) 
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Particularly, they showed that the zeroes of F{x) (for G = 0) play the roll of equilibrium positions for the soliton. In 
the case that only one zero, xq exists, the condition for stability for the kink/antikink is: 



dF 



dx 



J X — X(i 



{ 



> for the kink. 
< for the antikink. 



(3) 



For the stable case, the inhomogeneity can trap the soliton, but in the case that F{x) possesses more than one zero, 
the stability condition may become too complex. Here the extended character of the soliton arises, and interesting 
phenomena appear such as the interaction between the structure of the soliton (which is not a point-like particle) and 
the inhomogeneities, and between the shape modes of the soliton themselves. The external force F{x) can change the 
spectrum of small oscillations about the soliton and additional bounded states can exist. 

For the soliton of the 0* equation without perturbation, there are only two bounded states (besides the continuous 
spectrum): the translational mode and the shape mode. For specific values of the parameters that define the force 
F{x) used in Ref. ||^, not only an increase of the number of shape modes can exist, but in certain cases these can be 
unstable. Moreover the continuous spectrum can lose stability and the soliton becomes unstable against interaction 
with phonons. 

The above mentioned authors also considered the problem with the time dependent force G{x, t) and found that 
if G(x, t) has a spatial shape such that it coincides with one of the eigenfunctions of the stability operator of the 
soliton, then it is possible to get resonance if the frequency of the force also coincides with the resonant frequency of 
the considered mode. This means, for example, that energy can be given to the translational mode (or any of the 
shape modes) using a G{x,t) coupled to the translational mode (or any of the shape modes). 

When F{x) has three zeroes, this is equivalent to a double well potential (like the one found in the Duffing equation 
1^,^) chaotic motion of the soliton is possible applying an additional periodic force G{x,t) for a determined set of 
values of the parameters. 

In this paper we take into account the extended character of both the soliton and the impurity. We show that these 
considerations lead to the existence of a finite number of soliton internal modes that underlies a rich spatiotemporal 
dynamics. 

We introduce impurities of the N{x)4> type, where N{x) is a function with a bell shape. An impurity of this kind 
but using delta functions has been presented in Ref. . In our paper we consider a finite width impurity and a finite 
width kink and show the striking differences between this approach and the traditional one (structureless solitons and 
delta-function- like impurities) . 

We present a model for which the exact stationary soliton solution in the presence of inhomogeneities can be 
obtained and the stability problem can be solved exactly. To achieve this purpose we solve an inverse problem in 
order to have external perturbations which are generic and topologically equivalent to well-known bifurcation model 



We choose the "exact" solution such that the differential operator that appears in the stability problem is a Posch- 
Teller potential that can be solved exactly. Besides, the "generated" external force and the impurities have important 
physical properties. In particular, the inhomogeneous force F{x) is equivalent to the pitchfork bifurcation canonical 
form ||3^ and the iV(a;)0 impurity is topologically equivalent to the S{x)(j) type impurity that is very frequently used. 

Furthermore, we demonstrate the sensibility of the soliton internal dynamics to the inhomogeneity width even for 
the isolated impurity case. 

Our paper is organized as follows. In Sec. II we present a description of our model and we give specific physical 
interpretations of the equations under consideration. In Sec. Ill we study the equilibrium positions of the soliton 
and its stability. We analyze the interaction of the soliton with the inhomogeneity created by the interplay between a 
finite- width and 0-dependent impurity and the already studied inhomogeneities independent of (/) ||^ . We also consider 
the action of time dependent forces fitted to the shape of the translational mode of the soliton. In Sec. IV we describe 
the numerical simulations that confirm the theoretical results. We use the Karhunen-Loeve decomposition to relate 
the excitation of the shape modes spectrum of the soliton with the bifurcations. In Sec. V we present the interaction 
of the soliton with radiation modes. Finally, in Sec. VI we summarize and discuss our results and also present 
some concluding remarks. In the Appendices we outline the numerical method and present the Karhunen-Loeve 
decomposition. 



The topological solitons studied in the present paper possess important applications in condensed matter physics. 
For instance, in solid state physics, they describe domain walls in ferromagnets or ferroelectric materials, dislocations 




II. THE MODEL 
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in crystals, charge-density waves, interphase boundaries in metal alloys, fluxons in long Josephson junctions and 
Josephson transmission hues, etc. ||Tl| , ^ 

Although some of the above mentioned systems are described by the (/)^-model and others by the sine-Gordon 
equation (and these equations, in their unperturbed versions, present differences like the fact that the sine-Gordon 
equation is completely integrable whereas the <p'^-mode\ is not) the properties of the solitons supported by sine- 
Gordon and <p'^ equations are very similar. In fact, these equations are topologically equivalent and very often the 
result obtained for one of them can be applied to the other |^ . 

Here we consider the (j)'^ equation in the presence of inhomogeneities and damping: 

0.. - ^tt - 70* + ^ - 0') = -N{x)(l) - Fix), (4) 

where F{x) is a function with (at least) one zero and N{x) is a bell-shaped function that rapidly decays to zero for 
X — > ±oo. 

In ferroelectric materials (f> is the displacement of the ions from their equilibrium position in the lattice, ^ (0 — 4''^) 
is the force due to the anharmonic crystalline potential. Fix) is an applied electric field, and N{x) describes an 
impurity in one of the anharmonic oscillators of the lattice ||35| . In Josephson junctions, (f> is the phase difference of 
the superconducting electrons across the junction, F(x) is the external current, and N{x) can describe a microshort 
or a microresistor g . In a Josephson transmission line it is possible to apply nonuniformly distributed current sources 
{F{x)) and to create inhomogeneities of type N{x) using different electronic circuits in some specific elements of the 
chain |ll|,|6|. 

In the present paper the functions F{x) and N{x) will be defined as, 

F{x) = ]^A{A^ - 1) tanh(Ba;), (5) 



^ ' 2 cosh2(Bx) ^ ' 

The case F = const, has been studied in many papers (see e.g. [p2[). Here Eq. (^ represents an external field (or 
a source current in a Josephson junction) that is almost constant in most part of the chain but changes its sign in 
a; = (this is very important in order to have soliton pinning [^). Microshorts, microresistors and/or impurities in 
atomic chains |^ are usually described by Dirac's delta functions {S{x)) where the width of the impurity is neglected. 
The function N{x) is topologically equivalent to a 6{x) but it allows us to consider the influence of the width of the 
impurity. 



III. STABILITY ANALYSIS 

Let us consider the Eq. (||) and assume the existence of a static kink solution (j>k{x) that corresponds to a soliton 
placed in a stable equilibrium state created by the inhomogeneities F{x) and N{x). 
Analyzing the small amplitude oscillations around the kink solution (j)k{x), 

4>{x,t)=^k{x)+i^{x,t), (7) 

we get, for the function ip{x, t), the following equation: 

^xx - ^tt - li^t + \{l - i4>l + 2N[x))ij = 0. (8) 

Studying of the stability of the equilibrium solution (pkix) leads to the following eigenvalue problem (we have 
introduced i(j{x,t) — f{x) exp(A<) into Eq. (||)): 

-fxx + l{i<l>l-l-2N{x))f = rf, (9) 

where F = — — 7A. 

Let us now study some particular cases: If F{x) = 0, and the function N{x) is described by the expression (|^), 
then it can be shown that the exact solution for the kink in equilibrium at the position x = is. 
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(j)k{x) =ta,nh{Bx), (10) 
and that the discrete eigenvalue spectrum is described by the following formula, 

Tn = B^{A + 2An-n'^ -2); (11) 

where the parameter A is defined as, 

A(A + l) = ij + 2. (12) 

The integer part of A ([A]) defines the number of modes of the discrete spectrum. 

From Eq. ( pl| ) the stability condition for the translational mode (n — 0) can be obtained: 

< 1. (13) 

It is worth noticing that if the coefficient of N{x) is negative, the equilibrium position created by the impurity of 
the N{x)(j) type is stable for the soliton, and that the stability condition is independent of the polarity of the soliton. 
The opposite occurs when a inhomogeneity like F(x) is considered. 

Furthermore, it is necessary to point out the differences between the case in which the soliton equilibrium position 
is due to a zero of F{x), and the case in which the kink is trapped in the effective well created by the impurity N{x)(j). 
In the former the characteristic width and the number of internal modes that can be excited are smaller than those 
for the free kink, and even smaller than those for the kink in equilibrium in an unstable position. While, in the latter 
the characteristic width and the number of internal modes that can be excited are greater than those for the free kink 
and for the kink in unstable equilibrium in a repulsive impurity. This is due to the inverse proportionality between 
the width of the impurity and B (just like the width of the kink). The impurity is more stable as B diminishes. 

Finally, let us consider the case in which both types of inhomogeneities {F{x) and N{x)4>) are present in Eq. (|^). 
These functions are defined as in equations (^) and (^. 

For these functions the exact solution describing the static soliton can be written: 

(jjkix) = Ata,-Dh{Bx). (14) 
The spectral problem (Eq. ^) brings the following eigenvalues for the discrete spectrum: 

^-A2-i + s2(A + 2An-n2-2); (15) 



2 2 



here A is defined as. 



A 
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A(A + 1) = — + 2. (16) 

The stability condition for the translational mode is, 

165-* + 25^(5- 7^2) + (1 - yl2)2 <o. (17) 

When this condition is not fulfilled (the equilibrium position a; = is unstable) and > 1, then there will exist 
three equilibrium points for the soliton: two stable (at points x = xi > and x — X2 < 0) and one unstable at point 
a; = 0. This happens because for large values of |a;| the leading inhomogeneity is F(x), which is non-local and not 
zero at infinity. This inhomogeneity acts as a restoring force that pushes the soliton towards the point a; = 0. As a 
result of the competition between the local instability induced by N{x)(j) at x = and the non-local inhomogeneity 
F{x), an effective double- well potential is created. This is equivalent to a pitchfork bifurcation. 

Also note that when ^ = then A = 1, i.e., there is an unstable mode created by the impurity N(x)(f) even for a 
flat initial condition. This contrasts with the A = result expected for a delta-like impurity and with the result of 
Gonzalez and Holyst |6j for the other kind of inhomogeneity {F{x)). 



IV. QUASIPERIODIC AND CHAOTIC SOLITONS 

In this section we present numerical results for Eq. (^ perturbed with inhomogeneous external forces, impurities 
and time-periodic forces. First, we show, for inhomogeneities of the F(x) type, the bifurcations leading to the chaotic 
regime. Later, we study quasiperiodic regimes as the soliton internal modes are excited by an impurity of the type 
N(x)(f). Finally, we present bifurcations as the internal modes are excited for the case in which the two types of 
inhomogeneities {F{x) and N{x)(t)) are present. 
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A. Scenario of inhomogeneous external forces 



We simulate the following equation, 

^xx - <l>tt - l(l)t + I {(l> - 0^) = -F{x) - G{x, t), (18) 
where static and time-dependent forces are inhomogeneous and given by the following expressions: 

F{x) = ]^A{A^ - l)tanh(Bj;) + ^A{4:B^ - A^) smh{Bx) cosh-^{Bx) (19) 

and 

G{x,t) = .cos(a;t) + __^^_L__^) . (20) 

The force F{x), firstly introduced in Ref. can be obtained as a result of an inverse problem for which 
an effective potential equivalent to the canonical form of the pitchfork bifurcation is desired. For A^ > 1 and 

B^ 1 + ~ ^ '^^ presence of a double-well potential for the soliton, (i.e., the model is a Duffing-like 

psf soliton oscillator). Note that a force with three zeroes does not imply necessarily the existence of a double- well 
potential for the soliton. 

We will attend to the dynamics of the center of mass as defined by Eq. (AS). Figure |l| presents the bifurcation 



diagram (A = y |, i? = y A = 2, w = 1.22, E = 0.1, xi = 2.5 and 7 = 0.450) as the amplitude j/ of the time- 
dependent driving force is increased. The period doubling cascade for low values of the control parameter corresponds 
to oscillations of the soliton in one well. At 1^ = 0.312 an unusual discontinuous transition interrupts the cascade: the 
system switches from a period-eight solution to a period-two solution but in the other well as can be appreciated in the 
time series presented in figures |2](a)-(b). At = 0.318 the period-eight solution is recovered and a two-band regime 
with n-periodic (Fig. ^(a)), quasiperiodic (as the two-torus attractor presented in Fig. ||(b)) and chaotic attractors. 
The quasiperiodic attractor is a spatiotemporal effect: it is generated by the activation of a shape mode of the soliton 
which provides a frequency that is incommensurate to the driving frequency. Figures ^ and ^(c) reveals for v — 0.332 
the onset of jumps of the soliton between the two wells. This regime corresponds to Duffing-like chaos for the soliton. 
Figure ^(b) also presents the Poincare map which reveals the high-dimensional chaotic motion of the soliton which 
can be ascribed to an increased activation of shape modes (which is not possible for a Duffing-like chaotic particle). 
The intense activity around Xc.m. = (Fig. |l|) is also due to the extended character of the soliton. We have verified 
that at v ^ 0.356 the soliton prefers its deformation instead of its destruction and its dynamics returns to a periodic 
solution (figures |l| and ||(d)). 

Figures0(a)-(b) present time series for the center of mass of chaotic solitons. In the first case (^ = 0.3 and 7 = 0.505, 
the rest parameters are preserved) the soliton jumps between two wells. When the damping is increased (Fig. ^(b), 
7 — 0.550), the soliton is constrained to move mainly in one well. Corresponding Poincare maps (Fig. ||(a)-(b)) show 
the contrast of dimensionality (and consequently, of the number of modes and/or of the effective number of degrees 
of freedom) of these chaotic motions. Figures ||(c)-(d) present the spatiotemporal evolution for these chaotic solitons. 
Figure ||(c) evidences the jumps of the domain wall between the two wells, whereas Fig. ||(d) shows such a domain 
wall oscillating in one of the wells. Notwithstanding, note also in Fig. ^(d) the appearance and disappearance of 
deformations (depicted by yellow zones) in the other well due to tunnelling of mass/energy. 



B. Scenario of finite-width impurity 

We simulate the following model, 

<l>xx - ^tt - 7'/'* + ^ - = -N{x)cb - G{x, t), (21) 
where external forces are time-dependent and inhomogeneous. Right-hand functions are given by. 
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and 



V cosjajt) 

One may think that sohton dynamics in presence of an attractive impurity must be simple. Notwithstanding, that 



is not the case: as has been discussed in Section [II when the width of the impurity is finite, a large number of 
internal modes can be excited. It can exist energy exchange between these modes and the translational mode bringing 
a complex dynamics. 

Figure ^(a) shows the phase space for a period-three solution {A = 1.0, B ~ 0.25, A — 3.772, uj = 1.0, ly — 0.16 
and 7 = 0.1). As the amplitude v of the time-dependent driving force is increased, the excitation of an internal mode 
can provide a frequency incommensurate to the driving frequency, generating a two-torus quasiperiodic attractor 
(Fig. 11(b) for — 0.195, the rest parameters are the same as the previous case). On the other hand, decreasing of 
the damping parameter 7 can provide even more incommensurate frequencies to the driving frequency as reveals the 
three-torus quasiperiodic motion presented in figures ||(c)-(d) for 7 = 0.01 and v = 0.12. 

C. Scenario of finite-width impurity and inhomogeneous forces 

We will consider the following model: 

<^.. - 0tt -l^'t + li^- 0') = -Nix)^ - Fix) - G(x, t), (24) 
where N{x) and F{x) are given by Eq. (|) and Eq. (|), and G{x,t) is given by the following expressions. 



G(x,t) = lycos(Lut) [ — —rr-, tt H nrrrr, tt • (25) 

^ '\cosh^{B{x-xi)) cosh^{B{x + xi))J ^ ' 

We will study the sequence of bifurcations as the driving amplitude v is increased and other parameters remain 
fixed [A = 1.22, B = 0.32, A = 2, w = 1.22, xi = 2.5 and 7 = 0.3). Figure ^presents the time-averaged spatial profile 
for different values of the amplitude v of the time-dependent driving force: — 0.20 (period one), ly = 0.28. ..0.45 
(quasiperiodicity) , v — 0.55 (chaos) and v — 0.60 (period one). 

Figures ||(a)-(d) present a sequence of "quasiperiodic bifurcations": torus entangles as the amplitude of the time- 
dependent driving force increases and shape modes are activated. 

Figures ^( a)- (b) present an chaotic soliton for v = 0.55. The Poincare map (Fig. ^(a)) reveals the high-dimensional 
chaotic motion due to the activation of many internal modes whereas Fig. |9|(a) presents the temporal evolution of 
the soliton and shows that the kink profile is still sustained. For these parameters the soliton is at the edge of its 
destruction due to the activation of the shape modes. We can consider this regime as fully developed shape chaos for 
the soliton. 

Along this paper we had emphasized the importance of the internal modes (or shape modes) of the soliton that can 
be excited as the soliton interacts with inhomogeneities. Therefore we perform a Karhunen-Loeve decomposition for 
the sequence of dynamic attractors already presented for the scenario of finite impurity and homogeneous force. 

Figure |l^(a) reveals the increasing excitation of the discrete internal modes as the system evolves into a chaotic 
regime as well as the sudden change of the spectra for the final state that correspond to solution in which periodic 
motion is regained {i/ — 0.60). The periodic solution for ly = 0.60 corresponds to the highest deformity of the kink 
profile (Fig. 0). This agrees with the higher contribution to the dynamics of the first two modes whereas all the 
rest of the modes decreased their contribution. Furthermore, the first shape mode replaces the translational mode as 
the leading mode of the dynamics. Figure p^(b) presents the leading Karhunen-Loeve eigenmodes for the period-one 
solutions that initiates and ends the sequence of bifurcations considered in this section. The eigenmode for v — 0.20 
appears to be the superposition of a pair of translational modes centered at the equilibrium points for the soliton. 
Similar situation occurs for v = 0.60 but the eigenvalue appears to be the superposition of a pair of shape modes. 
Figures |ll|(a)-(b) the striking difference of the temporal evolution of the period-one solitons for v — 0.20 and v — 0.60. 

V. INTERACTION OF THE SOLITON WITH RADIATIVE MODES 

Any process that involves inelastic interactions or accelerations of the soliton leads to the emission of quasi-linear 
waves (radiation). This phenomenon occurs by means of the modes of the continuous spectrum (radiation modes [p9|) 
that for the case of the unperturbed equation are given by 
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/fc(a;) = e''=^ atanh^ (I) -6ifctanh(|) - (l + 4/c2) . (26) 

The interaction of a soliton with an inhomogeneity results in an emission of radiation that can be calculated using 
the method of McLaughlin and Scott ||2^ that relies on the construction of a Green's function that consists of a 
bilinear combination of eigenfunctions (i.e., in our case, of the eigenfunctions that we have already presented). The 
radiation problem also has been addressed using other perturbative methods [|| . 

In this paper we have focused our attention on the stability of the translational mode. Note that when the soliton 
is in an equilibrium position created by the impurities that is stable for the translational mode, the soliton as a whole 
is stable against the emission or absortion of radiation (these modes are usually called phonon modes). Moreover, in 
this case the kink is oscillating in an effective potential well for which radiative effects are exponentially small 

Notwithstanding, we remark that when the equilibrium position is unstable the shape modes can also become 
unstable (under certain conditions) and this can destroy the soliton. Moreover, under a certain condition the modes 
of the continuous spectrum can also be unstable. This is very surprising as it leads the soliton to be unstable against 
the emission and the absortion of radiation modes. This interesting effect is shown in Figure |l2| where we present the 
interaction of a radiation mode with the soliton under instability conditions for the continuous spectrum. The usual 
methods for calculating the radiation can not give this result. For instance, the phenomenon shown in Figure |l2| is 
produced by the evolution of the system 

(t>xx - <Ptt ~l<l>t + \{c^- 0^) = -F{x), (27) 

F{x) = 5(45^ - 1) tanh(Sa;) (28) 
under initial conditions that represent the supperposition of a radiation mode and a soliton, 

0(x, 0) = 2B ia.-ah{Bx) + C {3 sva{Bx) tanh(Ba:) - 6 cos{Bx) tanh(Ba;) - 5 sm{Bx)} , (29) 



Mx,0)=0 (30) 

For AB^ — 1 < the translational mode is unstable. For IQB^ — 1 < the first shape mode also becomes unstable. 
Moreover, for 125^ — 1 < the soliton becomes unstable against the interaction with radiation modes (Fig. [l^ . 

In Figure we present the case AB^ > 1. Here the amplitude of radiation mode is greater than the case presented 
in Fig. Notwithstanding, the original shape of the kink is recovered. 

The continuous spectrum (which concerns with radiation) was considered in our calculations and numerical exper- 
iments. In the numerical experiments there are emission of radiation that is absorbed by the resistance (dissipative 
terms). 

When the soliton is trapped the translational mode has a relevant role in contrast with the case for which the free 
soliton enters in interaction with an impurity. For the last case the shape modes indeed play a determinant mode. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have presented the importance of considering the soliton as an extended particle as it interacts 
with real inhomogeneous media. 

We have shown that a finite-width impurity can activate a large number of soliton internal modes. In fact, for an 
impurity of the N{x)(j) type with a stable point the number of mode increases (as the width of the impurity increases) 
in comparison with the free soliton case. Surprisingly, this contrasts with the inhomogeneous external force F{x) with 
a stable point case, where the opposite occurs. An interesting finite-size effect occurs when both inhomogeneities, 
F{x) and N{x)(j), exhibiting only one equilibrium point when considered individually, generate a double-well effective 
potencial for the soliton as they interact between each other and with the soliton. 

In addition, we have predicted surprising effects when more than one impurity is considered. For instance, the 
stability condition depends on a variety of parameters that we can control (i.e., the width and the height of the 
impurity as well as the width of the soliton) and which we present as very relevant to the whole dynamics. 

We have shown the importance of the internal modes of the soliton as they can generate shape chaos for the soliton 
as well as cases in which the first shape mode leads the dynamics. 

Our approach has been also applied to the stochastic resonance of solitons that can oscillate in a double-well 
potential p^ . Our approach provides the possibility to encounter resonances linked to the different internal modes 
of the soliton as well as to tune the time-dependent driving force to a selected mode. 
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The results presented in this paper are very relevant for the concrete physical systems presented in Sec. II. as 
they concern the pinning problem (i.e., the stability problem of the soliton in an equilibrium point created by the 
inhomogeneities) . As we have seen, it is not trivial to determine whether the equilibrium position is stable or not 
for the soliton when we are in the presence of external fields and inhomogeneities having finite widths and, also the 



extended character of soliton is considered. Real experiments ^9 40 1 where this kind of phenomena has been observed 
show that, in certain cases, the description of the impurities using delta functions can lead to erroneous conclusions. 

Soliton oscillators in Joscphson junctions have been studied intensively in the last years due to their application 
as sources of radiation ||4l[| . In this context the study presented in this paper of the sustained oscillations of a soliton 
in an effective potential well has great importance. Gr0nbech- Jensen and Blackburn have studied a system of 
coupled Josephson junctions as a super-radiant power source. Their point is that the oscillator ceases to be a point-like 
oscillator and exhibits an spatial degree of freedom. This leads to a radiation power greater than the expected from 
the point-like oscillators theory. In this paper we have shown that under the effects of the inhomogeneities the soliton 
exhibit an internal spatial degree of freedom. Furthermore, there is an increased number of activated internal modes. 
Particularly, we have shown that for impurities of the type N{x)(j) with a stable equilibrium position, the hyperradiant 
effect could be achieved varying the impurity width. In this case the width of soliton will be greater than the width 
of the unperturbed soliton. This is also valid for impurities of the type N{x) sin which can describe a microshort in 
a sine-Gordon-like systems. The chaotic oscillations and the soliton explosion phenomanon (due to soliton instability 
against the interaction with radiation modes) can be also important for the design of soliton based devices. 
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APPENDIX A: NUMERICAL METHOD 

We have integrated our t/f^-like equations using a standard implicit finite difference method with open boundary 
conditions 03,(0, t) = (j)x{l,t) = and a system length I = 80. 
We use a kink-soliton as initial condition: 



{x, 0) = Atanh 



B{x ~ xq) 



(Al) 



and 



0t(x,O) 



-ABvn 



cosh ■ 



B{x - Xq) 



(A2) 



here A and B are constants whereas Xq = 1.0 and t;o = are the initial position and velocity of the center of mass of 
the soliton. 

We define the position of the center of mass of the kink-soliton as. 



1-1% "^"t^ldx 



(A3) 



APPENDIX B: THE KARHUNEN-LOEVE DECOMPOSITION 



The Karhunen-Loeve decomposition |^-B^ allows to describe the dynamics in terms of an adequate basis of 
orthonormal functions or modes. The field w(x, t) to be decomposed represents the fluctuations of 0(a;, t) with respect 
to the time-averaged spatial pattern (j)ave{x). We find a basis of orthonormal functions 5'„(x) by solving an integral 
equation whose kernel is the two points correlation function K{x, x') = {uj{x, t)u!{x' , t)) (here (...) means time average). 
The functions 5'„(a;) are the eigenfunctions of the integral equation. 
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/ K{x,x')^r^{x')dx' = Xr,^r.{x). (Bl) 

Jo 

The eigenvalues A„ can be regarded as the weight of the mode n. 
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FIG. 1. Scenario of inhomogeneous external forces. Bifurcation diagram for the position of the center of mass of the kink 
sohton {u = 0.26.. .0.36, 7 = 0.450). 



FIG. 2. Time series Xc.m. vs t for the position of the center of mass of the kink soliton, 7 = 0.450. (a) Period eight, 1/ — 0.310. 
(b) Period four, v = 0.312. (c) Chaos, u = 0.328. (d) Period three, u = 0.360. 

FIG. 3. Periodic, quasiperiodic and chaotic attractors, 7 = 0.450. (a) Phase space Vc.m. vs Xc.m., v = 0.322. (b) Poincare 
map reveahng a two-torus (y = 0.327) and a strange attractor (y — 0.332). The inset presents the detail of one of the toruses. 



FIG. 4. Time series Xc.m. vs t for the position of the center of mass of the kink soliton, 1/ — 0.300. (a) Chaotic jumps of the 
soliton between the two wells, 7 = 0.505. (b) Chaotic soliton motion mainly in one well, 7 = 0.550. 

FIG. 5. Chaotic solitons, u = 0.300. (a) Poincare map revealing a high dimensional attractor as the soliton jumps between 
the two wells, 7 — 0.505. (b) Low dimensional attractor for kink motion mainly in one well, 7 = 0.550. (c)-(d) Spatiotemporal 
evolution for previous attractors. 

FIG. 6. Scenario of finite- width impurity, (a) Phase space Vc.m. vs Xc.m. for a period-three solution, 1/ — 0.160 and 7 = 0.100. 
(b) One-torus quasiperiodic solution, u — 0.195 and 7 = 0.100. (c) Phase space for a three-torus quasiperiodic solution, 
f = 0.120 and 7 = 0.010. (d) Poincare map for the three-torus attractor. 

FIG. 7. Scenario of finite-width impurity and inhomogeneous forces, {0(a;)) vs x for 1/ = 0.20 (period one), u = 0.28. ..0.45 
(quasiperiodicity) , 1/ = 0.55 (chaos) and 1/ = 0.60 (period one). 7 = 0.300. 

FIG. 8. Scenario of finite-width impurity and inhomogeneous forces: quasiperiodic bifurcations, 7 = 0.300. (a)-(b) Phase 
space and Poincare map for 1/ = 0.280. (c)-(d) Phase space and Poincare map for f = 0.350. 



FIG. 9. Shape chaos, 7 = 0.300 and f = 0.550. (a) High dimensional Poincare map. (b) Spatiotemporal evolution. 

FIG. 10. (a) Karhunen-Loeve spectra for the sequence of bifurcations presented in Fig. |^. (b) Plot of the first mode of the 
Karhunen-Loeve spectrum for 1/ = 0.20 and v — 0.60. The inset shows the traslational mode and the first shape mode. 

FIG. 11. (a)-(b) Spatiotemporal evolution of the period one solutions that initiates and ends the sequence of bifurcations 
presented in Fig. Q. 

FIG. 12. (a)-(j) Interaction of the soliton with radiation modes: soliton explosion. After the explosion an anti-kink is formed. 
Note that the initial kink profile is stretched because the equilibrium position is unstable whereas the opposite occurs for the 
final stage. (L = 80, 7 = 0.1, B = 0.2 and C = -0.001). 

FIG. 13. (a)-(b) Interaction of the soliton with radiation modes under stability conditions for the continuous spectrum. The 
amplitude of the radiation mode is greater than the case presented in Fig. |l^ as is evident in the snapshot (a) for t = {L = 80, 
7 = 0.1, B = 0.8 and C = -0.011). 
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